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ABSTRACT 



Context. Difference imaging has proven to be a powerful technique for detecting and monitoring the variability of unresolved stellar 
sources in M 31. Using this technique in surveys of galaxies outside the Local Group could have many interesting applications. 
Aims. The goal of this paper is to test difference imaging photometry on Centaurus A, the nearest giant elliptical galaxy, at a distance 
of 4 Mpc. 

Methods. We obtained deep photometric data with the Wide Field Imager at the ESO/MPG 2.2m at La Silla spread over almost 
two months. Applying the difference imaging photometry package DIFIMPHOT, we produced high-quality difference images and 
detected variable sources. The sensitivity of the current observational setup was determined through artificial residual tests. 
Results. In the resulting high-quality difference images, we detect 271 variable stars. We find a difference flux detection limit corre- 
sponding to m R ^24.5. Based on a simple model of the halo of Centaurus A, we estimate that a ground-based microlensing survey 
would detect in the order of 4 microlensing events per year due to lenses in the halo. 

Conclusions. Difference imaging photometry works very well at the distance of Centaurus A and promises to be a useful tool for 
detecting and studying variable stars in galaxies outside the local group. For microlensing surveys, a higher sensitivity is needed than 
achieved here, which would be possible with a large ground-based telescope or space observatory with wide-field imaging capabilities. 

Key words. Galaxies: individual: Centaurus A - Galaxies: stellar content - Stars: variables: general - Gravitational lensing 



1. Introduction 



Using difference imagin g techniques (Tom anev & Crotts| [l996; 
Gouldl 1 996t I Alardl2000h it is possible to detect and monitor stel- 
lar variability in highly crowded fields where the individual stars 
are unresolved. During the past decade these techniques have 
been used to this aim by se veral microlensing surveys towards 
the M agellanic Clouds (e.g. lAlcock et alj|200fj ; [Tisserand et all 
2007 ) and the Andromeda galaxy (M 31) (e . g. iRiffeser et alj 
20031 ICalchi Novati et alj 120051: Ide Jong et al.ll2006l) . This has 
led to the successful detection of several microlensing events and 
tens of thousands of previously unknown variable stars in M3 1 . 
Following these results the question arises whether the applica- 
tion of these techniques can be extended to even more distant 
objects. Both the detection of variable stars and of microlensing 
in galaxies outside the local group would have several interesting 
applications. 

The number of known variable stars in galaxies could be 
multiplied, which would be interesting in itself, and also useful 
for, for example, distance determination and the study of stellar 
populations and star formation histories. Microlensing surveys 
outside the Local Group would also be an important tool to study 
the luminosity functions and the halo compact object content of 
galaxies other than the Milky Way and M 3 1 . On the one hand 
this would be useful as a further measure of the contribution of 
massive compact halo objects (MACHOs) to galaxy halos. On 
the other hand such microlensing surveys could be an important 
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step to the generalisation of the picture obtained from the sur- 
veys in the Local Group. Both the Milky Way and M31 are two 
large spiral galaxies in a relatively low-density environment. It 
would be interesting to also study the compact object content 
in the halos of galaxies in other environments or in the halos of 
other types of galaxies. 

If the dark matter does indeed have a baryonic MACHO con- 
tribution, this component could be different in elliptical galax- 
ies compared to spirals, since they are believed to have under- 
gone different formation histories and star formation episodes. 
Furthermore, elliptical galaxies tend to live in more dense envi- 
ronments, where galactic encounters are more common and old 
stellar populations would have been pushed out to larger radii 
into the halo. Microlensing surveys in elliptical galaxies would 
also enable the study of these stellar halos, that might be impor- 
tant repositories of faint stellar mass. 

A first attempt at a microlen s ing su rvey in an elliptical 
galaxy was made by iBaltz et ail d2004l) . Using observations 
taken with the Wide Field and Planetary Camera 2 (WFPC2) 
on the Hubble Space Telescope, they performed a search for mi- 
crolensing events in the giant elliptical gala xy M 87. Due to the 
large distance to M 87 of about 16 Mpc dTonrvet al.ll2001l) do- 
ing this ex periment proved difficult in practice. In their study, 
IBaltz et al.1 d2004 find only seven variable sources, 1 of which is 
consistent with microlensing, although it is bluer than expected 
for a typical microlensing event. 

In this paper we present a pilot study of a ground-based vari- 
ability study using difference imaging photometry towards the 
giant elliptical galaxy Centaurus A (Cen A, NGC 5128). Apart 
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Table 1. Summary of data used 



Date 


Exposures 


Total t exp 


Seeing 






(sec) 


(") 


April 17 2005 


4 


4000 


0.99±0.06 


April 30 2005 


9 


9000 


0.83±0.10 


May 9 2005 


8 


8 000 


0.78±0.03 


May 11 2005 


9 


9000 


0.84±0.13 


May 12 2005 


9 


9000 


0.92±0.07 


May 29 2005 


9 


9000 


1.00±0.14 


May 30 2005 


9 


9000 


0.82±0.10 


May 31 2005 


9 


9000 


0.80±0.13 


June 1 2005 


9 


9000 


1.02±0.15 


June 8 2005 


9 


9000 


0.79±0.07 



from the general interest of a variable star and microlensing 
study in an elliptical galaxy discussed above, this survey would 
have additional scientific interests. From the warped gas layer, 
the outer isophotes and the kinematics of the planetary nebu- 
lae, there are indications that t he potential and therefore the halo 
around Centaurus A is triaxial dHui et al.ll 19951) . This means that 
the lines-of-sight toward different parts of the galaxy have dif- 
ferent path-lengths through the halo. Because of this, the mi- 
crolensing event rate will be asymmetric over the face of Cen A, 
much like the microlensing rate due to halo lensing is asymmet- 
ric between the near and far side of the disk of M 3 1 . The spatial 
distribution of halo microlensing events could therefore also be 
used to constrain the shape o f the dark halo of Cen A. Because 
the distance to Cen A (4 Mpc. lTonrv et alJuOOll) is much smaller 
than the distance to M 87, the sensitivity will be much higher and 
the crowding of variable sources much less strong. Nevertheless, 
Cen A lies five times further away than M 31, so the project is 
not without challenge. 

Section [2] describes the data and methods that were used in 
this study. In Section [3~2l the resulting difference images are de- 
scribed and the detected variable sources presented. We deter- 
mine the detection efficiency for variable sources in Section [33l 
The prospects of ground-based (halo) microlensing surveys to- 
wards Centaurus A are estimated in Section [4] based on our re- 
sults and a simple model of the halo of Centaurus A. Finally, in 
Secfion[5]we present our conclusions. 



2. Data and methods 

Goal of the observational set-up was to obtain deep, 9 000s, pho- 
tometry in the broad R c filter during each of a series of nights 
spread over a several month period, resulting in both high sensi- 
tivity and sufficient temporal coverage to detect faint variable 
sources. Observations were done during 12 nights spread be- 
tween April 15th and Jun e 8th 2005 with the Wide Field Imager 
(WFI; lBaade etalJll999t) at the ESO/MPG 2.2m telescope at 
La Silla, Chile. With its field-of-view of 34'x33'(pixel scale 
0.238"per pixel) the WFI enables the monitoring of Cen A in 
one single pointing. The total exposure time per night was di- 
vided over 1 000s exposures to help with cosmic ray rejection. 
A small subset of data was discarded because of poor seeing 
(>1.3"), leaving data for 10 nights, and 9 nights with the total 
exposure time of 8 000s or longer. A summary of the exposure 
times and average seeings is listed in Table[TJ The seeings listed 
are the average of all exposures from the four central chips in the 
array, with the quoted errors the standard deviation. 



Table 2. Mean zero points and their uncertainties 



Chip 


Zero Point 


Uncertainty 




(mag) 


(mag) 


ccd50 


24.238 


0.040 


ccd51 


24.334 


0.048 


ccd52 


24.113 


0.040 


ccd53 


24.312 


0.035 


ccd54 


24.311 


0.029 


ccd55 


24.363 


0.028 


ccd56 


24.293 


0.034 


ccd57 


24.171 


0.029 



The data reduction was performed with the Astronomical 
Wide-Field Imagin g System for Europ43 (AstroWISE, 
IValentiin et al.ll2007h . A short description of the main steps of 
the data reduction process is given here. 

- Bias subtraction: 

For each night a master bias image was created from be- 
tween 10 and 42 raw bias frames after doing an overscan 
correction and outlier rejection. These master bias frames 
were subtracted from all overscan-corrected science frames 
to remove 2-D bias patterns. 

- Flat-fielding: 

Master dome and twilight flats were constructed for each 
night from between 5 to 20 exposures. The master dome and 
twilight flat-fields are combined into a master flat-field used 
to correct the science frames. 

- Astrometry: 

The astrometry was derived using routines from the Leiden 
Data Analysis Center (LDAC). A preliminary astrometry is 
first estimated as a simple shift of the 8 chips together to 
align the field with the coordinates in the USNO A2.0 cata- 
logue. The final astrometry for each chip is derived by using 
2nd order polynomials to fit the position of stars detected in 
all exposures to their USNO coordinates. 

- Photometry: 

Since the photometry of the images is scaled to a reference 
image during the difference imaging procedure, only mean 
zero points were computed in AstroWISE to convert the 
ADU flux of variable stars to magnitudes. The zero points 
were derived as the median of the zero points obtained from 
the Landolt standard fields observed during the whole pe- 
riod of the observations assuming a mean extinction in R c 
of 0.09 mag. On average, 4 standard fields were observed 
per night, with an average of 100 standard stars landing on 
each chip. The zero points and their uncertainties are listed in 
Table [2] The uncertainties are the standard deviations of the 
zero points derived from all 40 standard field observations. 

- Resampling: 

Finally, images were regridded to a pixel scale equal to 
0.2"pix _1 using the LANCZOS3 algorithm in SWARlQ 

2. 1 . Difference imaging photometry 

For the difference image analys i s the same pipeline is used as 
was employed by Ide Jong et al.l d2006l) in their analysis of M 
31 data. This method is based on the procedures for matching 

1 A detailed description of the system can be found at the AstroWISE 
portal: www.astro-wise.org/portal 

2 For information on SWARP see http://terapix.iap.fr/ 
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Fig. 1. Example of an area in Centaurus A and its difference image. Left: detail of the reference image showing the dust lanes and 
star forming regions to the south-east of the centre of Centaurus A. Right: detail of the difference image of June 8th 2005 showing 
the same region of Centaurus A. 



the point spread function (PS F) between two images described 
by [Tomaney & Crottsl 0996) and although it has a different 
approach to the problem of PSF-match ing of two images from 
the method developed by lAlardl d2000t) it produces difference 
images of similar quality. The difference image pipeline runs 
in IRAF and makes use of the DIFIMPHOT package written 
primarily by Austin Tomaney. Taking the reduced images as 
input, the pipeline goes through the following steps. 



- Variable sources are detected on the difference images 
using SExtractor dBertin & Arnoutsll 19961) . where we select 
groups of 3 adjacent pixels that are at least 3<x above or 
below the local background. Bright stars, saturation spikes 
and bad pixels are masked out during this step. 

- Using PSF-fitting photometry the brightness of all variable 
sources is measured from the difference images and light- 
curves are constructed. 



- The images are aligned to a common astrometric reference 
frame to an rms precision better than 0. 1 pixels using stars 
spread over the whole image. 

- Using cuts on the seeing, astrometric shift and sky bright- 
ness, the best ~15 exposures are selected and combined into 
a high signal-to-noise reference image. 

- All images obtained during the same night are combined 
into one image. These nightly image stacks will provide 
the photometric measurements for the detection of variable 
sources and their light-curves. Thus, since we use data 
obtained during 10 nights we have 10 temporal epochs in 
the subsequent analysis. 

- After photometric calibration and matching of the PSF 
the reference image is subtracted from each nightly stack, 
resulting in a difference image for each night. The PSF of 
the nightly stacks and the reference image is measured from 
bright but unsaturated stars on the image. From these PSFs 
a convolution kernel is derived in Fourier space to convert 
the better seeing image to the seeing of the worst seeing 
image before subtraction. To cope with PSF variations 
across the chips, each chip is divided in 16 subregions, for 
each of which the PSF matching is done separately, using 
a minimum of 10 stars. Photometric calibration is done 
using a subset of non-variable bright stars. The resulting 
difference images are dominated by shot noise in which 
variable sources show up as positive or negative residuals. 



3. Results 

3. 1 . Difference image quality 

The difference images that result from the above described pro- 
cedure are of very high quality, with noise levels close to the 
theoretical photon noise. In Figure Q] a cutout from the refer- 
ence image and the corresponding cutout from the difference 
image of June 8th 2005 are shown. While most of the differ- 
ence image shows only noise, variable stars in the star forming 
regions leave positive and negative residuals. Bright foreground 
stars also leave residuals in the difference image due to imperfect 

subtraction. 

As shown theoretically by iGouldl dl996l) . it is in principle, 
possible to create difference images where the noise is domi- 
nated by photon noise. Noise due to systematic effects such as 
PSF variation and geometric and photometric alignment can, 
when handled carefully, be brought down to levels where they 
cannot compete with the photon noise. To assess how well our 
difference images approach the limit of pure photon noise, we 
have determined the noise in several regions of highly vary- 
ing background intensity where no obvious variable sources are 
present. In the upper panel of Figure [2] these noise levels are 
plotted versus the background intensity for four different images 
with different FWHM. There is a strong correllation with higher 
noise in brighter regions, which is a clear signature of the pho- 
ton noise. Apart from this, there are offsets between the differ- 
ent epochs, which depend on a combined accuracy of the PSF- 
matching, photometric and astrometric alignment, seeing varia- 
tions between individual exposures etc. 
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Fig. 2. Difference image noise levels for a set of regions with 
different background intensity. Top: standard deviation of pixel 
values for four epochs with different seeing values: May 11 
2005 (0.84", crosses), May 12 2005 (0.92", triangles), June 1 
2005 (1.02", squares), and June 8 2005 (0.79", circles). Bottom: 
Same as in top panel, but scaled to the theoretical photon noise. 
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Fig. 3. Residual strength, in terms of theoretical photon noise, of 
a set of bright, resolved stars of different magnitudes for four 
observation epochs. Top: Residuals in a high surface bright- 
ness region (p,R ^28 mag/arcmin 2 ) in the difference images of 
May 11 2005 (short dashed), May 12 2005 (long dashed), June 
1 2005 (dot-dashed), and June 8 2005 (dotted). Bottom: Same 
as top panel, but for a low surface brightness region {/ur ^29 
mag/arcmin 2 ). 



The lower panel of Figure [2] shows the same, but now di- 
vided by the theoretical photon noise. For a single exposure, the 
photon noise is simply the square root of the number of counts 
in a pixel; in this case both the reference image and the nightly 
images are stacks of several exposures, so that in each the pho- 
ton noise should be the square root of the counts divided by the 
square root of the number of exposures going into the stack. The 
quality of our difference images turns out to be very good, with 
the noise very close to the theoretical photon noise. There might 
be a slight degradation at very high background levels, because 
very close to the centre of Centaurus A the number of stars that 
can be used to construct an accurate PSF model is low, due to 
saturation of the chips. 

Another interesting test is to look at the residuals around 
resolved, bright, but unsaturated stars. Because of their higher 
brightness, secondary effects such as imperfections in the PSF- 
matching kernel or photometric scaling will show up more 
strongly. Figure [3] shows the noise divided by the theoretical 
photon noise for stars of varying brightness in a high back- 
ground region (pg ~ 28.0 mag/arcmin 2 , upper panel) and in a 
low background region (pn ~ 29.0 mag/arcmin 2 , lower panel). 
In both panels the residuals are much higher than the photon 
noise for bright stars, and decreasing to the photon noise level 
at faint magnitudes. In the high background region this decrease 
is faster because the photon noise is higher. That the residuals at 
the bright end are higher in the high background region is caused 
by the poorer PSF modelling near the centre of Centaurus A, due 
to the lack of sufficiently bright, unsaturated stars. It should be 
noted that this problem does not affect image subtraction meth- 
od s based on direct pixel-t o-pi xel fits, such a s the one developed 
bv lAlard & Luptonl £1998) and lAlardl (120001) . 



3.2. Variable sources 

Residuals were detected in the difference images using 
Sextractor and at all positions with at least two detections light- 
curves were built. A by-eye inspection of light-curves and dif- 
ference images to reject spurious detections due to bad pixels or 
foreground stars, resulted in a catalogue of 271 variable sources. 
A catalogue with the positions of the variables is published in 
electronic form with this article. 

The majority of these variables are expected to be Cepheids 
and long period variables (LPVs) such as Mira variables, since 
these classes are both common and bright. There is also one 
Nova in the sample. That we would catch one Nova during its 
peak is not unexpected considering the Nov a rate of ~28 per 
year for Centaurus A (Ciar dullo et alJlT990t) . With ten points, 
the sampling of the light-curves is too sparse to allow an accu- 
rate period and amplitude determination for the Cepheids. For 
the LPVs the situation is worse since their periods are longer 
than the total time span of our observations. We did, however, 
apply two different period finding algorithms to the light-curves, 
namely the multi-harmonic p eriodogram method described by 
ISchwarzenberg-Czernvl(11996|) and the normalised periodogram 
method developed bv lLombl (1 1976 ). both of which are designed 
to work on unevenly sampled data. For 29 variables the periods 
determined with the two methods lie within 10% of each other 
and between 2 and 50 days. These limits exclude periods that are 
longer than the observing period and the aliases that occur close 
to one day since this is the minimum spacing between the data 
points. These reasonably robust periods, which are in the range 
of typical Cepheid periods, are listed in the electronic catalogue 
of variables. Of the other 242 variables, the majority of the ones 
in the star forming regions seem chaotic. These are most likely 
periodic variables for which the periods cannot be determined 
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Fig. 4. Selection of light-curves from variable stars in Centaurus A, showing difference flux in ADU/s as a function of the time in 
days from the start of the observations on April 17th 2005. (a) Bright Cepheid-like variable with a period of 43 days, (b) Same 
variable as in (a) but showing two periods of the folded light-curve, (c) Cepheid-like variable with a period of 2.8 days, (d) Same 
variable as in (c) but showing two periods of the folded light-curve, (e) Nova variable, erupting on May 29 and reaching its maximum 
around May 31 2005. (f) Long-period variable reaching its minimum at the beginning of June 2005, with the maximum occurring 
outside our time range. 



because of poor phase coverage. The majority of the variables 
away from the centre are monotonically increasing or decreas- 
ing and can be classified as LPVs. 

In Figure |5]the positions of all 271 variable sources are plot- 
ted. There is a clear concentration of variables in the star form- 
ing regions in the dust ring in the central parts of Centaurus 
A. The majority of the stars with confirmed periods between 
2 and 50 days are located in these regions. This distribution is 
not unexpected, since Cepheids are massive stars with masses 
of at least 3 or 4 times that of the Sun and up to several tens 
of solar masses, and therefore quite young. On the other hand, 
LPVs belong to old stellar populations and are expected to be 
more evenly distributed. Also shown in Figure [5] ar e the loca- 
tions o f the fields monitor e d for variability bv lReikuba et al.l 
d2003l) and iFerrarese et ail d2007l). Us ing ISAAC at the Very 
Large Telescope (VLT), Reikuba et al.l d2003l) were able to de- 
tect 1504 LPVs in an area just over 10 square arcminutes. None 
of these are detected in our survey, due to the lower sensitiv- 
ity of our data and the fact that LPVs are sign ificantly fainter 
at opti cal wavelengths compared to the infrared. Ferrarese et al.l 
d2007l) used the Wide Field and Planetary Camera 2 (WFPC2) 
on board the Hubble Space Telescope to detect Cepheid vari- 
ables, resulting in a sample of 56 Cepheids and 70 other periodic 
variables. Due to the difference in sensitivity and the fact that 
Ferra rese et al.l (120071) use clear periodicity as one of their detec- 
tion criteria, we only identify three stars that are both in our and 
in their sample. For completeness, also the Nova variables from 



the extragalactic part of the General Catalogue of Variable Stars 
(lArtvukhina et al.lll996l) are plotted. To give the reader an im- 
pression of the quality of the light-curves, we present the light- 
curves of two periodic variables, a Nova and an LPV in Figure 
|4] The Nova variable shown in panel (e) of the figure has a peak 
difference flux corresponding to to«=20.2, making this by far 
the brightest variable in our sample. A period-amplitude plot is 
shown in Figure[6]for the 29 variables with robust periods. The 
amplitudes were determined by taking the difference between 
the highest and the lowest point of the lightcurve; these am- 
plitudes estimates are lower limits since the actual maxima and 
minima might not be sampled by the data. The dashed line in- 
dicates the period-amplitude relation found for type I Cepheids 
in M 3 1 using the same techniques as employed here (Ide Jong! 
120051) . 



3.3. Detection efficiency 

To determine the detection efficiency for variable sources we put 
artificial residuals in the difference images and try to "detect" 
them. Of course, there are different ways to detect sources in 
variability studies such as pixel lensing surveys. The method 
used in this paper to detect variable stars searches the differ- 
ence images for groups of pixels that exceed a certain thresh- 
old above or below the background. This method was also used 
by the MEGA collab oration in their M 31 microlensing sur- 
vey (cf . Ide Jong et al. 2006). In another frequently used method, 
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Fig. 5. Positions of the 271 variable sources detected in our sur- 
vey, with respect to the centre of Centaurus A. Closed circles 
indicate variable sources for which a reasonably accurate period 
could be determined . The two squares outline the fields in which 
Reik uba et ail d2003l) detected a total of 1504 long -period vari- 
ables using ISAAC at the VLT; the WFPC2 field-of-view used by 
iFerrarese et alJ (120071) in their survey for Cepheids is located just 
to the right from the centre of the Fig ure. Crosses indicate th e 
positions of Nova variables taken from Artvukhi na et al.1 (fl996). 



light-curves are built at each pixel position in the image, after 
whic h the light-curves are ide ntified that show a variable signal 
(cf. [Calchi No vati et"aTll2005l) . To be as general as possible we 
will use both approaches in our detection efficiency analysis and 
show that the difference in sensitivity between the two is small. 

We select a subset of six difference images for the detection 
efficiency tests that are of high quality and from CCDs in the 
centre of the array. The images are from the nights of April 30, 
May 30 and June 8 2005 and correspond to chips 2 and 7 from 
the CCD array. This way reductions in sensitivity due to poor 
seeing and PSF degradation in the perifery of the focal plane are 
avoided. Although these kind of effects are unavoidable in real 
surveys, the detection efficiencies we retrieve will represent up- 
per limits for the current observational setup. Because the back- 
ground from unresolved stars varies strongly in intensity across 
the field, the detection efficiency will vary spatially. To probe 
this effect we divide both chips in the part where the background 
surface brightness is lower than /i«=28.5 mag/arcmin 2 and the 
part where it is higher, and consider these regions separately. 
The high background region includes the small but very bright 
centre of Centaurus A. We randomly place residuals of a fixed 
(difference) flux in the selected difference images, 1 000 at a time 
in the low background regions and 250 at a time in the smaller 
high background regions. Fluxes probed are 2, 1.5, 1, 0.75, 0.5, 
and 0.25 ADU/s, corresponding to 23.6, 24.0, 24.4, 24.7, 25.2 
and 25.9 magnitudes; in total we put 20000 residuals with these 
fluxes in each region, respectively. The actual PSFs measured 
from the stacked nightly images, which are divided in subre- 
gions to deal with possible PSF variations across the field, are 
used for building the residuals. Photon noise is included in the 
residuals. 
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Fig. 6. Periods and amplitude estimates of the 29 variables with 
robust period determinations. The dashed line represents the ap- 
proximat e period ampli tude relation found for type I Cepheids 
in M 31 dde Jongll2005l) . scaled to the distance of Centaurus A 
and corrected for the different gains. 



After inserting the residuals, SExtractor is run on the differ- 
ence images and groups of 3 adjacent pixels that are at least 3<x 
above the local background are selected. Comparison of the list 
of detected sources with the input list of residuals gives the de- 
tection efficiency of this residual detection method. At each posi- 
tion where an artificial residual is inserted, we also perform PSF- 
fitting photometry to measure the flux at that position. We con- 
sider the residual to be detected if the recovered flux is at least 
3 cr above zero, where the cr is the photometry error. This error 
is calculated by the PSF-fitting photometry routine and based on 
the statistics of the fit and includes a correction for photometric 
scaling inaccuracies during the image subtraction. The fraction 
of sources detected with this PSF-fitting photometry gives a sec- 
ond measurement of the detection efficiency. 

Figure [7] shows the detection efficiencies for both detection 
methods and all regions. There is indeed a difference in the de- 
tection efficiency in the bright regions (upper panels) and the 
fainter regions (lower panels), although not very large. In both 
cases the detection efficiency is 100% down to a flux of 1.5 
ADU/s and virtually 0% at 0.5 ADU/s, with the transition more 
gradual in the bright regions. This more gradual slope is due to 
the larger variation in the background intensity. These results 
show that even very close to the bright centre of Centaurus A, 
residuals with a flux of 1 .5 - 1 ADU/s are still detected with high 
efficiency. There is not much difference between the SExtractor- 
based detection method (left panels) and the photometry-based 
method (right panels). A variability survey with the WFI at the 
ESO/MPG 2.2m telescope with 9 000s exposure time per epoch 
thus will have a detection limit of at best slightly below 1 ADU/s. 
This corresponds to ~24.5 trig, or for a variable of absolute mag- 
nitude Mr--\ to a magnitude difference in the order of 1 magni- 
tude. Since the majority of Cepheids has amplitudes of less than 
1 magnitude this means that only the very bright ones or ones 
with high amplitudes are detected. 
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Fig. 7. Detection efficiencies in percent as function of difference 
flux in ADU/s. Dashed lines show the results for chip 7 in the 
WFI array, dotted lines those for chip 2, where each line corre- 
sponds to a different epoch. The left panels show the efficiencies 
for the case where residuals are detected as groups of significant 
pixels in the difference images, the right panels for the case of 
PSF-fitting photometry. In the upper panels the results for the 
regions with a background surface brightness higher than 28.5 
mag/arcmin 2 are plotted, in the lower panels for the regions with 
lower surface brightness. 



4. Implications for microlensing surveys 

Now that we have shown that Difference Image Photometry 
works well at the distance of Centaurus A, we want to investigate 
the prospects of ground-based microlensing surveys towards this 
galaxy. In the previous section the detection efficiency for a sin- 
gle residual of a certain flux in a high quality difference image 
was determined. Taking this as the detection efficiency for a mi- 
crolensing event of the same maximum amplitude is of course 
a gross oversimplification. To be able to "detect" a microlens- 
ing event, enough information must be obtained from the light- 
curve to be reasonably certain of the lensing nature of the de- 
tected variability. In practice one would want to detect the event 
in more than one difference image and with sufficient signal-to- 
noise to recover certain characteristic signatures, such as light- 
curve shape and achromaticity. Even assuming a survey with ob- 
servations every night, a detection limit of 1 ADU/s is optimistic, 
but can serve as an upper limit. 

In the light of the arguments discussed in Section[T]we con- 
centrate on microlensing by objects in the halo of Centaurus A. 
If we assume the halo is a cored isothermal sphere with velocity 
dispersion <x and core radius r c , the density is given by, 



a~ 

p{r) = — 



1 



2nG r 2 + r 2 ' 



(1) 



By definition, a source is considered to be microlensed if it is 
lying within the Einstein radius of a lens, i.e. if it is magnified at 
least a factor of 1 .34. The microlensing optical depth r is defined 
as the number of lenses within one Einstein radius of the line-of- 



sight towards the source and therefore gives the fraction of stars 
that is microlensed at a given time. This means that t is given by, 



■r 

Jo 



ni(r)nR\ dr — 



f b p(r) 4nGM leas D LS D 0L 



Dos 



dr, (2) 



where «; is the number density of lenses, fb is the fraction of 
halo mass in the lenses, Re is the Einstein radius, and Dol, Dos 
and Dls are the distances between observer, lens and source. We 
integrate from the position of the source to the outer radius of 
the halo r/,. For external galaxies like Centaurus A, Dol — Dos 
and Z\s — r, so that, 



r(r ) 



2f b cr 2 



r 2 + r 2 l + r 2 



dr, 



(3) 



for a source star located at a projected distance ro from the centre 
of the galaxy. Integrating gives, 



T(r ) 



fbo 1 



In 



r 2 + r 2 



(4) 



For a MACHO mas s fraction ft, of 1 0%, a halo velocity disper- 
sion of -100 km/s dPeng et al.ll2004T) and reasonable values for 
r/,, ro and r c this gives values in the order of 10~ 7 , meaning that 
at any time one in ten million stars is being microlensed. 

However, at the distance of Centaurus A we are in the regime 
where only very high amplification events ar e detected, the so- 
called spike pixel-lensing regime defined by iGouldl d 19961) . In 
this regime the impact parameter f3, expressed in units of the 
Einstein Radius Re, is very small, ft <k 1, and inversely pro- 
portional to the maximum amplification A: f3 ~ 1 1 A. Suppose a 
source has to change in brightness by a certain amount Ld et in or- 
der to be detected. This means that for a microlensing event to be 
detected, the maximum impact parameter is /3 ~ 1 /A - L/L d et 
for a source star with luminosity L. The number of detectable mi- 
crolensing events at any time is now given by the optical depth 
multiplied by the square of the impact parameter, 



N det ^ t ■ f3 2 (L) = t J" dLf(L) f2 



L 2 
L 



(5) 



det 



where f(L) is the luminosity function of the source stars. Using 
the fact that the luminosity of the galaxy L ga i = J dLf(L)L 
and the fluctuation brightness of the stellar population L = 
J dLf(L)L 2 1 L ga i we obtain, 



Mlet - I'- 



LL, 



gal 



(6) 



The absolute magnitude o f Centaurus A is Mr =s -22 
dLauberts & Valentiinlll989h and a typical value for the fluctu- 
ation brightne ss of an old stellar population is in the order of 
R - -1.3 (e.g.|BuzzQnfl993; To nry et al.l l2001). Given our de- 
tection limit of niR ^24.5, corresponding to Mr ^-3.5 this leads 
to a number of detectable microlensing events at any time of 
A^det = 10 6 5 t. Assuming a MACHO halo mass fraction of about 
10% this means that at any moment ~0.3 microlensing events 
can be seen. 

The rate F at which microlensing events due to halo lens- 
ing occur is given by the number of events A^et divided by the 
average duration of the events, measured by the time during 
which the magnification is above half the maximum magnifi- 
cation ?fwhm- In the regime where 1, we can approximate 
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^fwhivi = V3y6?E (iGondololl 19991) . Here fg is the Einstein cross- 
ing time, the time it takes the lens to cross the full Einstein disk, 
given by, 



t E = 



2Re 

v± 

,1/2 



2 
v± 



4GM lens D LS D 0L 



Dos 



i8xl0 v(*)" ! (-^)" ! pi 

\M Q ) \10kpcj \l00km/s 



(7) 



(8) 



where v± is the transverse velocity of the lens with respect to the 
observer-source line of sight. A typical value of ?e for 1 M halo 
objects in Centaurus A is therefore in the order of a year. The 
event rate is now given by, 

r - r f dLf{L)^ Itvmm - V3- f dLf(L)-^- 



V3^ 

?E ^det 



(9) 



, and filling in the numbers gives r a a/3 ■ 10 74 T/f E . A gain 
assuming an ft, of 10% we get an expected detectable event rate 
of ~4 per year. The average event duration is simply N^/T ~ 1 
month. 

An added complication for studying microlensing due to 
MACHOs in the halo of Centaurus A is that these microlensing 
events have to be separated from different possible contaminants, 
such as variable stars and lensing by other lens populations. Most 
periodic variables can be recognised as such by extending the 
time baseline of the survey sufficiently. Cataclysmic variables 
are not period, but can be identified based on the lightcurve 
shape and colour evolution. Consi dering the expected Nova rate 
for Centaurus A of ~28 per year (ICiardullo et al.l l 1990). which is 
several factors higher than the halo microlensing rate, the need 
for sufficiently dense time sampling and multi-band photometry 
is clear. Microlensing by stars in Centaurus A, rather than by 
objects in the halo, will also have to be disentangled from the 
halo lensing signal. As for M 3 1 and M 87, this so-called "self- 
lensing" will domi nate over halo lensing i n the central parts o f 
Centaurus A (cf. iBaltz et al.ll200l 120041: iRiffeser et ai] 120061) . 
The cleanest halo lensing sample is therefore obtained in the 
outer parts of the galaxy. 



5. Conclusions 

In this paper we have demonstrated that difference imaging 
techniques can be successfully applied to ground-based data of 
Centaurus A, a giant elliptical galaxy at a distance of 4 Mpc. 
Using the Wide Field Imager at the ESO/MPG 2.2m telescope 
at La Silla, and a total exposure time of 9 000s per epoch, we 
reach a difference flux detection limit corresponding to a bright- 
ness of Mr ^24.5. We have detected and photometered 271 vari- 
able sources in Centaurus A, mostly Cepheids and long-period 
variables and a Nova. The majority of these variables were pre- 
viously unknown and are available as an online catalogue with 
this paper. Clearly, these techniques can be used as powerful 
tools to enlarge the samples of known variable stars in and to 
study stellar populations of Centaurus A and other galaxies in 
the nearby universe. At these depths the difference images do 
not suffer from crowding, and can be photometered straightfor- 
wardly. 

The microlensing optical depth towards Centaurus A due to 
lenses in the halo is of the order of t ~ 10~ 7 if we assume that 
10% of the halo mass is in compact objects. Given this optical 



depth and our detection limit the expected microlensing event 
rate is ~4 per year over the whole galaxy. At large distances it 
will be impossible to identify the source stars of individual mi- 
crolensing events and blending by stars and variables that are 
too faint to be detected will play a more important role. Also, 
only very high amplification events can be observed, m eaning 
that t he microlensing is completely in the spike regime dGouldl 
1996). Even more than for microlensing in M 31, the information 
that can be obtained for individual events is limited. Statistical 
analysis of event rates, spatial distribution and timescale distri- 
bution is the only way of extracting scientific information out of 
these data. This again emphasises the need for large samples of 
microlensing events. 

The pilot experiment described here shows that it would 
be possible to put useful limits on the number of microlensing 
events in Centaurus A, but a survey that is several magnitudes 
deeper would be required. This could be achieved with a large 
ground-based telescope equipped with a wide-field imager large 
enough to take in all of the galaxy in one exposure, or better 
still, with a wide-field s pace observatory. We therefore endorse 
the plea of lBaltzl d2005l) for inclusion of microlensing of nearby 
galaxies into the observing programmes of upcoming wide-field 
cosmological survey facilities. 
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